Aging of the frictional properties induced by temperature variations 
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The dry frictional contact between two solid surfaces is well-known to obey Coulomb friction 
laws. In particular, the static friction force resisting the relative lateral (tangential) motion of solid 
surfaces, initially at rest, is known to be proportional to the normal force and independent of the area 
of the macroscopic surfaces in contact. Experimentally, the static friction force has been observed 
to slightly depend on time. Such an aging phenomenon has been accounted for either by the creep 
of the material or by the condensation of water bridges at the microscopic contacts points. Studying 
a toy-model, we show that the small uncontrolled temperature changes of the system can also lead 
to a significant increase of the static friction force. 

PACS numbers: 46.55.+d Tribology and mechanical contacts, 65.40.De Thermal expansion; thermomechan- 
ical effects, 62.20.Hg Creep. 



I. INTRODUCTION. 

Granular materials, or more generally macroseopie 
solids in frictional contact, at mechanical equilibrium are 
a priori considered as athermal systems, meaning that 
thermal agitation has no significant feedback on the me- 
chanical degrees of freedom at the grain scale. However, 
uncontrolled temperature variations inevitably produce 
dilations of materials and one can wonder whether they 
can alter the properties, in particular the mechanical sta- 
bility, of such systems. 

Several recent studies showed that temperature cycles, 
even of small amplitude, can induce the slow compaction 
of dry granular materials [ll-Q- The mechanisms, in par- 
ticular the role played by the confining walls, are still 
under debate but one can anyway conclude that the tem- 
perature variations cause aging: the properties of the ma- 
terial, here at least the density, evolve in an irreversible 
manner with time. 

In granular matter, several sources of aging were iden- 
tified. First, for "dry" granular materials in a humid at- 
mosphere, the condensation of microscopic liquid bridges 
at the contacts between grains, leads to a cohesive force 
Q : the nucleation of the bridges being an activated pro- 
cess, the mechanism leads to a logarithmic increase of the 
angle of avalanche Q and of the static frictional co- 
efficient 0, Q with time. Second, for immersed granular 
materials, chemical reactions at the surface of the grains 
results in soldering them to each other Q. Again, the 
angle of avalanche @ and the static frictional coefficient 
0,131 observed to increase logarithmically with time. 

In the solid friction [HI, which characterizes the 
mechanical contact between solid surfaces, a similar ag- 
ing phenomenon also takes place: the static frictional 
coefficient is also observed to increase logarithmically 
with time for contacts between various materials like 
metals [lH, rocks [l^, Bristol paper Again, for 

"dry" friction in a humid atmosphere, the condensation of 
liquid bridges at the micro-contacts between the flat, but 
nevertheless rough surfaces, can account for a part of the 
phenomenon |15| . However, when effects of the humidity 



are suppressed, aging is still observed. The phenomenon 
is due to the creep of the material at the micro-contacts, 
as proven by the dependence of the aging dynamics on 
the overall temperature of the system [l6| . Experiments 
even exhibited a logarithmic aging of the micro-contacts 
themselves [l7| . 

To our knowledge, the possibility for the temperature 
variations to be responsible for a part of the aging ob- 
served in solid friction was still not evaluated. We can 
guess that the dilation or contraction of the materials 
which result from a change in the temperature alter 
the stress distribution between the micro-contacts and 
that, as a result, the system can evolve with time due 
to a pinning-depinning dynamics. Such dynamics of the 
micro-contacts has been proposed to explain the dissipa- 
tion associated with the relative motion of the surfaces 
and, thus, to account for the dynamic friction [l8l - [20j but 
effects of dilations have not been evaluated. The study of 
the hysteresis cycle of the tangential force induced by a 
quasi-static cyclic displacement is a related but slightly 
different problem 21|: small displacements, even small 
enough not to be considered as macroscopic sliding, in- 
duce irreversible micro-slips and a slow evolution of the 
system with time . The temperature variations would 
have a slightly different effect as they rather induce ran- 
dom stress variations, correlated at long range. 

In the present article, we propose the study of a mini- 
mal model mimicking the frictional contact between two, 
nominally flat, solid surfaces, a slider on top of a horizon- 
tal substrate. We first consider that, due to roughness, 
the real contact between the surfaces reduces to a large, 
but finite, number of microcontacts, themselves belong- 
ing to a small number of mesoscopic, coherent, contact 
regions [1^. Such coherent regions are composed of a 
large number of microcontacts, so that their contact with 
the substrate obeys the Amonton-Coulomb law for static 
friction. However, the number of microcontacts in a co- 
herent region is finite 



241 so that the associated static 



friction coefficient, /ig, is sensitive to the local microstruc- 
tures in regard and, thus, depends on the position on the 
substrate. In addition, coherent regions are elastically 



2 



connected one to another by the material the shdcr and 
the substrate are made of. 

In order to mimic the practical situation in the sim- 
plest model, we further reduce the problem to the study 
of a system consisting of two sliders connected by a linear 
spring. Each of the sliders accounts for one of the co- 
herent regions in frictional contact with a flat substrate 
whereas the spring accounts for the elastic coupling be- 
tween them. We shall demonstrate that, when the system 
is subjected to thermal dilations, the mechanical stability 
of the slider, accounted for by an effective static frictional 
coefficient, exhibits a significant increase with time due 
to the distributed values of the static friction coefHcient 



II. THE MODEL 

A. Description 

The system (mass M), subjected to its own weight Mg, 
lies on a flat and horizontal surface. In a simple approach, 
we further assume that each contact region sustains half 
the weight of the system, so that the normal force act- 
ing on one slider is = Mg/2. We introduce the rest 
distance between the two contact regions at the tem- 
perature Tq. Due to the elasticity of the slider material 
(assumed to be much softer than the substrate) , the small 
sliders are connected by a spring of length Zq and stiffness 
k and, being identical, they are associated with the same 
inertia, i.e. the same mass m = M/2. With these as- 
sumptions, the system reduces to two sliders, connected 
by a spring, moving on a flat surface as sketched in the 
figure [TJ Wc denote by xi and X2 the positions of the 
small sliders I and 2 respectively such that the length of 
the spring I = X2 ~ xi. 



coefficient, /Zd, which quantifies the rate of energy dissi- 
pation. We further assume, in agreement with standard 
observations, that /Xd < /^s- 

The dynamics of the system is induced by temperature 
changes which, due the thermal dilation of the materials, 
lead to changes in the natural length of the spring I ac- 
cording to: 

l{T)^lo[l + K{T-n)\ (1) 

where T stands for the temperature and k for the thermal 
expansion coefficient of the slider material. We assume, 
for the sake of simplicity, that the substrate does not 
dilate, so that the positions xi and X2 are not altered by 
the temperature changes, and that the spring constant 
does not depend on the temperature. Thus, due to the 
thermal dilation, the spring expands or contracts so that 
the slider 2 is subjected to the tangential force 

Fi^2 = -k{x2 -xi^l) (2) 

where we assume the x-axis to be oriented from the slider 
1 to the slider 2. In the same way, the slider 1 is sub- 
jected to F2^i = —Fi^2- If l-P'i-i-21 exceeds one of the 
static friction forces, fis.iN or fis.2N^ the corresponding 
slider starts moving and the system rearranges. After the 
slider has stopped, the corresponding value of the static 
frictional coefficient has changed because the slider lies 
at a different position on the substrate. The new value of 
the frictional coefficient is drawn at random from a distri- 
bution p(/is)- In addition, we assume for consistency that 
fig is always larger than the dynamic frictional coefficient 
/^d, namely p(^s) = for < Md- 

The aim of this study is to account for the evolution 
of the system properties, thus of fis,i and /is. 2, as time 
elapses (aging). 

B. System of equations 
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FIG. 1. Sketch of the considered situation. 

At rest, the contact between each slider and the sub- 
strate is characterized by the static frictional coefficient 
/is, which accounts for the value that the horizontal force 
must overcome to move the slider. Due to the heterogene- 
ity of the surfaces in regard, the local static frictional 
coefficient might take, at random, different values, /is,i 
and /is. 2, for the two contact regions. By contrast, we 
characterize the dynamical friction by a single constant 
value of [lA. Indeed, the surface explored by the sliders 
being homogeneous in average, the dissipative interac- 
tion with the substrate can reasonably be characterized, 
over the whole surface, by a single dynamical frictional 



The problem reduces to the study of two sliders, con- 
nected by a spring, in frictional contact with a substrate. 

First, one can write the equations governing the dy- 
namics of the sliders when in motion. Introducing the 
thermal dilation Q = k{T — Tq), one gets 

771X1 = -~k[xi - X2 + lo{l + e)] - fidNS{xi) (3) 

777X2 = —k[x2 — Xi - lo{l + 9)] — ^dNS{x2) 

where S denotes the sign function. The dynamics is char- 
acterized by a typical time scale Tdyn = \/rn/k. We re- 
mark that, in practice, Tdyn ~ lo/cs where Cg stands of 
the speed of sound in the material the macroscopic solid 
is made of. For a typical size Iq ~ 10 cm and usual values 
of Cs (about a few kilometers per second), we estimate 
Tdyn ~ 10-6 - 10-5 s. 

Second, the static solid-friction is accounted for by the 
coefficients figj (j = 1,2) such that, initially at rest, the 
slider j starts moving and, thus, the system rearranges 
when 

k\x2-xi-loil + 0)\^ fisjN. (4) 
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C. Rearrangements 

From now on, we assume that after n rearrangements, 
the sHders are at the positions and X2 associated with 
the values /x";^ and //"a of the static frictional coeffi- 
cients. From the condition one notes that the re- 
arrangements involve the weakest slider, namely the one 
associated to the smallest value of the static frictional 
coefficient. For convenience, we denote: 

Mmin = minl/^M ' <2) = (5) 

where we introduce the index z„ of the corresponding 
slider. With these definitions, we get the value 0" of the 
dilation 6 at which the (n + 1)*^ rearrangement occurs: 

_ n 

9:^-^^~l + fi-^,,NSie) (6) 

'0 

where the dimensionless normal force TV = N/{klo) has 
been introduced to lighten the notations. Eq. ^ holds 
true when the system dilates {6 > 0) or contracts {9 < 0), 
the sign S{9) being evaluated before the rearrangement, 
for e 9^. 

For 6 = 0", the slider i„ moves and reaches a novel 
static position such that: 

A.T,;„ = 2(-l)^" - fid)mSi9). (7) 

One can easily show that the second slider necessarily 
remains at rest so that the rearrangement induces a dis- 
placement: 

Axg = (-1)'" - ^ld)NloSi9) (8) 

of the center of mass G of the two sliders. 

In conclusion, when the dilation 9 reaches the critical 
value 9" , the weakest slider moves by Axi^ and reaches 
a novel static position whereas the other slider remains 
at rest. As a consequence, the static frictional coefficient 
fis,in takes a new random value, fJ,"f^, associated with 
the new position x"^^ = xf^ + Axi^^ . The position of 
the other slider and the corresponding static frictional 
coefficient remain unchanged. 

III. NUMERICAL ANALYSIS 

We shall consider the temporal evolution of the system 
when subjected to aleatory temperature changes. In the 
section IIII Al we describe how the temperature changes 
are accounted for and, then, in the section llll BI we report 
the behavior of the system. For practical purposes, we 
consider in this section a Gaussian distribution p{ps), 
namely 

p(Ms) = ^^exp (9) 



Although this distribution does not strictly speaking sat- 
isfy the condition p{fis) = for fis < A*d, this condition 
will be satisfied in practice as we restrict our study to 
the case cr^ ^ (/Jg — yUd)- 

A. Temperature changes 

We shall assume that the temperature T of the system 
fluctuates around the temperature Tq such that 9 is dis- 
tributed according to a Gaussian distribution "0(0) with 
zero mean and variance ag. In addition, due to its ther- 
mal inertia, the system exhibits a thermal characteristic 
time Tth which limits the dynamics of the temperature 
changes and, thus, of the dilations. We assume that the 
temperature variations occur on time scales much larger 
than the typical time scale of the rearrangement dynam- 
ics, that is, Tth 3> Tdyn. 

In order to account for the finite value of Tth, we con- 
sider the evolution of the system in the following way: 
At time tq = grth {q integer), the system is described 
by the positions and X2 , the frictional coefficients /.i" ^ 
and /i"2 a'^d the dilation 9{tq). A new value of 9{tq+i) is 
drawn at random from the distribution ip{9). Then, we 
consider ^ as a constant over the time interval [tq,tq+i], 
namely 

o='Jh±ii^m (10) 

Tth 

and one can calculate 0" from Eq. ([5]). The correspond- 
ing positions a;"^^ and Xj"*"^ after the rearrangement are 
calculated according to Eq. (O, and the smallest fric- 
tional coefficient among ^"^^ and ^"2 drawn anew 
from p(/is), the other remaining unchanged. The process 
is repeated while 9^' < 9{tq+i) if > 0, or 9'^ > 9{tq+i) 
if ^ < 0. When the evolution stops, the dilation between 
tq and tg+i has led to r rearrangements. The system has 
reached a new state and the process can be iterated, by 
choosing randomly a new value 9{tq^2)- 

B. Numerical results 

When the system is subjected to temperature changes, 
the values of the static frictional coefficients continuously 
change in time. As the rearrangements only involve the 
weakest slider, the maximum of the values /i"]^ and 

Mmax = max(/i^i,/i^2), (11) 

cannot decrease, and thus on average continuously in- 
creases with time t (Fig. [2]). We indeed observe that /imax 
increases rather logarithmically with time t and that, for 
instance, the aging is faster for larger ae ■ 

In addition, due to the random motion of the center 
of mass associated with the rearrangements (Eq. [7]), the 
system diffuses. We report the mean square displace- 
ment (xq) as function of time t (Fig. [3]) and observe that 
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FIG. 2. Maximum frictional coefficient /imax vs. dimensionless 
time t/rth - The is observed to increase almost 

logarithmically with time t, the aging being more effective 
when the amplitude ao of the temperature variations is larger 
(/Xd = 0.5, -p^ = 0.6, cr^ = O.OI, N = 10~'). 



log(t/rth), the diffusion being faster for larger ag. 



Thus, a dilation or a contraction of amplitude A9 leads 
to a typical number of rearrangements r « lAfll/l^"^^ — 
d"\. To average over the fluctuations of A9, wc define it 
quantitatively as the difference between two values and 
6' randomly drawn from ip{0). Hence A6 is a Gaussian 
random variable with zero mean and variance 2(Tg. We 
thus find for the average value (r) 



(13) 



Comparison with the results of numerical simulations 
(Fig. |4]) essentially confirms this prediction. A careful 
analysis of the numerical data however reveals the pres- 
ence of a small shift. The average number of rearrange- 
ments is found to be described by the phcnomenological 
form 



(j0 - ^idN 



(14) 



The small offset fi^N is intuitively expected to result 
from the absence of rearrangements when 9 changes sign. 
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FIG. 3. Normalized mean square displacement {x'q) / fj'gl'i 
vs. dimensionless time f/rth - The mean square displacement 
{xq) increases logarithmically with time t, the diffusion being 
faster when the amplitude ag of the temperature variations is 
larger (/id = 0.5, Jl^ = 0.6, cr^ = 0.01, N = 10"''). 
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FIG. 4. (a) Number of rearrangements n vs. time t - The 
number of rearrangements n increases, to within small fluc- 
tuations, linearly with the time t. (b) Average number of 
rearrangements < r > during Tth vs. amplitude erg - In aver- 
age, the slope dn/dt or equivalently < r > depends linearly 
on the amplitude ag according to Eq. (|14p . 



In the next section, we shall discuss further these re- 
sults theoretically, describing analytically how the aging 
and the diffusion depend on the parameters of the prob- 
lem at long times. 



From now on, we shall consider the evolution of the 
system in terms of the number of rearrangements n. The 
temporal evolution of the system can then be assessed by 
considering the simple relation between n and time t. 



IV. THEORETICAL ANALYSIS 

A. Dilation and rearrangements 

In a first step, we consider the number r of rearrange- 
ments that occur during the time interval rth. Consid- 
ering the Eqs. © and ([7]), the difference between two 
successive values of the dilation at which rearrangements 
occur can be estimated as 

\9:+'^9:\^2{jI^-f,,)NS{e). (12) 



B. Aging of the static frictional coefficient 

It is particularly interesting to focus on the dynam- 
ics of /iJJ,ax ^ this quantity plays an important role in 
the stability of the system (Sec. |V|. Actually, one sim- 
ply needs to notice that the dynamics consists in draw- 
ing randomly two values ^ and 2 the initial step 
n = I, and then drawing at each step a new value of 
the frictional coefficient to replace the smallest one. As 
a results, the value /x^j^^ is the maximum value of the 
set of n -|- 1 random values drawn independently. This 
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mechanism is actually very close to the one appearing in 
a standard model ofglassy dynamics, namely the Barrat- 
Mezard model [lH, [2^ . The statistics of /ijjjax then boils 
down to an extreme value problem of independent and 
identically distributed random variables, for which stan- 
dard results are well-known [l^,!!!]. In particular, the 
distribution of /^"^^x belongs for large ?? to one of the three 
classes of extreme value statistics, namely the Gumbel, 
WeibuU or Frechet class. However, we are here more in- 
terested in the typical value of ^^^x rather than by its 
relative fluctuations. This typical value can be estimated 
through the following simple scaling argument. We de- 
note as F{^s) the complementary cumulative probability 
distribution 

/•oo 

F{fi,) = / p{f,)dpi. (15) 

Having /^JJ^^x smaller than a given value fi is equivalent to 
the fact that the n + 1 random values drawn dynamically 
are below the value /i. The corresponding probability is 
simply [1 — i^(/i)]"+^. The typical value ^typ of /^^ax then 
satisfies 

[l-F(Mtyp)r+^«i (16) 

which for large n takes the form 

1 - F(/.typ) « « 1 - i (17) 

so that finally F{fityp) w 1/n. If the upper wing of the 
distribution p(/i) decays as p(/^) ^ e""^^ where c and a 
are positive constants (for instance, a = 2 for a Gaussian 
distribution), then i^((u) behaves in the same way as p{fJ-) 
up to algebraic prefactors. It follows that 

/ityp«p, + a^c-i/"(lnn)i/" (18) 

to leading order in n. Thus, in the case of the Gaussian 
distribution considered in the section [nil the simple the- 
oretical argument presented here predicts an increase of 
Mmax like 



/imax«7is+'7My21n(^-^^, (19) 

taking into account the correspondence n w {r)t/Tth be- 
tween the time t and the total number n of rearrange- 
ments (Eq. ITi)) . 

C. Slowing down diffusion 

In Sec. mil we obtained a logarithmic subdiffusion of 
the center of mass, which we now would like to interpret. 
We first observe that the evolution of the system can be 
divided in successive periods during which the strongest 
slider remains the same, that is, the frictional coefficient 



of the weakest slider is repeatedly drawn anew, without 
exceeding the value /^max of the strong slider. During 
such periods, the position of the weakest slider changes, 
due to the rearrangements, and the net displacement £ is 
limited by the extension of the spring. One can estimate 
£ ^ agio. When the chosen frictional coefficient becomes 
larger than the former value /imax (which can be called 
a record breaking), the weakest and strongest sliders get 
exchanged. Hence the mean square displacement of the 
center of mass after n rearrangements should be of the 
order of £^ times the number of record breakings in this 
sequence of rearrangements. 

It is a classical result from statistical record theory 
that the mean number of records that occur when succes- 
sively drawing n independent and identically distributed 
random values is, to leading order, equal to Inn for large 
n. Further using the relation n w {r)t/Tth, one eventu- 
ally finds for the mean square displacement of the center 
of mass 

{:cl)^<jlllln(^^y (20) 

thus accounting for the logarithmic time dependence and 
the order of magnitude of the numerical data reported in 
Fig. [31 

V. DISCUSSION 

We have shown that the dilations associated with the 
temperature changes are likely to induce the aging of 
the physical properties of the system (the maximum fric- 
tional coefficient increases with time and the diffusion 
slows down). At this point, it is particularly interesting 
to discuss, first, the orders of magnitude of the temper- 
ature changes necessary for the aging process to be at 
stake and, second, the possibility to assess the aging pro- 
cess experimentally. 

Let us now consider the order of magnitude of the di- 
lations likely to induce the aging. From equation (|14p . 
given that fid ^ 1, we deduce that the minimal ampli- 
tude of the dilations necessary to rearrange the system 
is of the order of N which can be estimated as follows. 
Consider that the two contact points, sketched by the 
slider 1 and 2, are located at the distance Iq one from 
the another. Taking into account the Young modulus, 
Y, the material is made of, one can estimate k ~ YA/Iq 
where A stands for the surface area of the slider cross 
section in the perpendicular direction. Considering the 
density p of the material and volume V AIq of the 
slider, one can estimate the normal force, due to gravity 
g, N ^ pAlog. Thus, we get N ^ pglo/Y. In practice, 
considering g = 10 m.s~^, p ~ 10^ kg.m~^, Y 100 GPa 
and the typical size = 10 cm, we get TV ~ 10^®. Thus, 
dilations of the order of 10~® are likely to rearrange the 
system. Such dilations, the typical thermal expansion 
coefficients being of about 10~^ K~^, correspond to tem- 
perature changes of about 10"'^ K. Thus, in usual cxper- 
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<r> t / T. 



FIG. 5. Variation (/^^ ~Ms)/''"m number (r)t/rth - We 
report data obtained for two different values of the dynamical 
frictional coefficient (Lines : = 0.5; Open squares : — 
0.4). The values of cr^ are chosen in accordance with the 
values of a^/{jl^ — ^^d) given in the figure (/J^ = 0.6, N = 
10-^ ag = lO-'^). 

imental conditions, the thermal dilations produce rear- 
rangements at the contact scale. Even more, it is not 
obvious to perform a control of the temperature that in- 
sures that rearrangements will not occur. 

Experimentally, one possible method to highlight the 
aging process consists in measuring the critical angle of 
stability: the substrate is slowly tilted until the system 
starts sliding as a whole for a critical tilt angle, ac (isj . 
In this configuration, the sliders are subjected to normal 
and tangential forces due to gravity and one can relate ac 
with an effective static frictional coefficient fi'f^ = tanac- 
We assessed numerically the temporal evolution of as 
described in the Appendix [A] We observe that the in- 
crease in yLtmax (Fig- 121) results in an almost logarithmic 
increase in fj.°^ with time t. As suggested by Eq. (|19p . we 
report the relative variation {fj,°^ — /Ig)/"'/^ ^ function 
of the typical number of rearrangements n w {r)t/Tth 
(Fig. [5]). We observe a collapse of the curves correspond- 
ing to the same value of cr^/ijls ~ Md), which suggests: 

^,fit)^Jl^ + aM^^^^). (21) 

From the numerical data, the function $ depends almost 
logarithmically on t at long times. Asymptotically, as a 
function of the number of rearrangements n(i), ^ 
/Ig -|- (/)(7^ (logn)'' with /3 of the order of unity and (j) a- 
constant which only depends on cr^/ ijl^ — Md)- The aging 
phenomenon is drastically enhanced when the width cr^ 
is increased. 

It is particularly interesting to estimate the expected 
order of magnitude of the effect. With TV ~ 10~® and 
typical values /Is ~ 0-6 and /id ^ 0.5 of the frictional 
coefficients, temperature changes of about 0.01 K corre- 
spond to ae ~ 10"'' and, thus, to typically < r >^ 50 
rearrangements during the characteristic time rth. From 
a typical thermal conductivity A of about 1 W/(m K) and 
a typical heat capacity C of about 10^ J/m^, one can esti- 



mate from the diffusion coefBcient A/C that rth ~ ^qC/X 
ranges from 10~^ to 1 s (indeed much larger than r^yn) 
for a typical size of the system ranging from 1 to 10 cm. 
Considering further a value cr^ ~ 0.02 (about 3 % per- 
cent of the average value Jig), we expect from Fig. [5] an 
increase of about 0.2 % of the effective frictional coeffi- 
cient after 10^ rth ~ 10^ — 10^ s. Note however that for 
CTp ~ 0.03 (about 5 % percent of the average value /Jg), 
the effect is drastically enhanced leading to a change of 
about 5 % after the same waiting time. Thus, small tem- 
perature variations are likely to induce, in a few hours, 
an increase of the effective frictional coefficient of about a 
few percents. Note that, from Eq. ((2T|) . we expect the ag- 
ing process to mainly depend on a^^ and ae (through (r)), 
the effect being expected to be larger for broader distribu- 
tion of the frictional coefficient and for larger amplitude 
of the temperature variations. Finally, it is important to 
notice that the thermal dilations are important because 
of the change klo9 they induce in the force between the 
two sliders and that the latter change is to be compared 
with the normal force N. Thus, a smaller k leads to a 
weaker effect of the temperature changes. This is prob- 
ably the reason why no or little aging is observed with 
soft materials [l5| . 



VI. CONCLUSION 

Reducing the study of the frictional contact between 
two solid surfaces to the study of two frictional sliders, 
connected by a spring, in contact with a flat substrate, we 
demonstrated that the uncontrolled thermal dilations of 
the system can lead to a significant increase in the effec- 
tive static frictional coefficient 11°^ with time. The evolu- 
tion of /imax at long times is assessed theoretically and the 
associated 11°^ is estimated numerically. Our simplistic 
model makes it possible to demonstrate that the process 
leads to an almost logarithmic increase in and we ac- 
count for the dependence on the parameters, in particular 
the distribution of the static frictional coefficient associ- 
ated with the solid-solid contact and the characteristics 
of the temperature changes (amplitude and characteris- 
tic time) . Considering typical orders of magnitude of the 
physical ingredients, we show that temperature changes 
can have significant effects, which can also contribute to 
the a ging of the frictional coefficient observed experimen- 
tally [ll]- 

Finally, one could wonder if the effect reported here 
would not be specific to the system we chose. We as- 
sumed that the contact between two flat and rough sur- 
faces reduces to a small number of coherent regions. On 
the one hand, one could claim that the contact between 
the surfaces rather consists in a large number of micro- 
scopic contact points and that Non-Amonton behavior 
of the microscopic contacts was revealed experimentally 
[30| . However, whatever the characteristics of the micro- 
scopic contacts (characterized by frictional mechanical 
properties or not), one can guess that the thermal dila- 
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tions, which induce a change in the distribution of the 
force network between the latter, lead to rearrangements 
of the less stable ones. The process again favors the most 
stable contacts and we expect the effective frictional co- 
efhcient to increase with time. On the other hand, even 
considering coherent regions exhibiting frictional proper- 
ties, one could claim that their number is likely to be 
significantly larger than considered in the present study. 
Again, we do not expect an increase in the number of 
contact regions to change qualitatively the aging process. 
However, the quantitative effect is difficult to predict the- 
oretically. The study of a 2D irregular network of sliders 
is thus pertinent and will be the subject of a forthcoming 
publication in which the aging, the diffusion, as well as 
the creep motion induced by a constant tangential force, 
will be considered. 



Appendix A: Effective static frictional coefficient 

In this appendix, we present the numerical method 
used to compute the time evolution of the effective fric- 
tional coefficient . To this aim, we consider that the 
substrate on which the two sliders lie can be tilted by an 
ajustablc angle a. 

The dynamics of the system on the horizontal sub- 
strate, as described in Sec. IIII Al is divided in time steps 
of duration Tth- Within each time step, the sliders evolve 
through a sequence of rearrangements during which the 
positions and static frictional coefficients change. Before 
each rearrangement, the values of /is,i, Ms. 2 and of the 
force -^1^2 ai'c copied to auxiliary variables, and the an- 
gle of avalanche corresponding to this precise set of vari- 
ables is determined through an iterative algorithm that 
we now describe. 

Due to gravity, when the substrate is tilted, the sliders 
on the incline are both subjected to the same positive 
force / = mg sin a. Correspondingly, the normal force 
now becomes = mg cos a. According to Eq. Q, the 
slider 1 (rcsp. 2) starts moving if one of the two following 
conditions are fulfilled: 

/ - Fi^2 > fis,iN, (Ala) 

/ + Fi^2 > t^s,2N. (Alb) 

One deduces from these relations that, when a is in- 
creased, the slider 1 moves first if 

2Fi^2 < (Ms.2-Ms,i)^, (A2) 

where Fi_j.2 is given by Eq. ([5]). The onset of motion 
occurs precisely when / = mg sin a reaches the value 
^ ^^ jN + Fi^2- If the condition ((M\ is not fulfilled, 
the slider 2 moves first, when / reaches the value /2 = 

Ats,2Af-Fi^2. 

The motion of the slider i induces the motion of the 
other slider if the condition ^ is satisfied, thus if (after 



simple algebra): 

|r >Md + ^(-l)'(Ms,i-A^s,2). (A3) 
Then, if / > HdN, the motion of the two sliders is accel- 
erated, and a macroscopic sliding of the whole system is 
observed. The corresponding value of tan ac = f / mg is 
recorded. In the opposite case, when / < /i^, the motion 
is damped, leading only to a short displacement of the 
system (creep). Once the system is again at rest, new 
values of fis,i and /Xs,2 are drawn from the Gaussian dis- 
tribution p(/Lts). The value of the force Fi^2 also needs to 
be recomputed, due to the change of the positions xi and 
X2. Since the determination of the new positions would 
require to integrate the equations of motion during the 
damped motion, which is time-consuming, we rather use 
a simple approximation. The new force Fi^2 is chosen at 
random, with a uniform probability, between the values 
/ — /is,i N and /is.2 N — f, where ^j.s,i and /is. 2 designate 
the new values of the static frictional coefficient. We are 
again led back to Eq. (|A1|) . and the process is then iter- 
ated until a macroscopic sliding is observed. 

If, however, Eq. (|A3p is not satisfied, the other slider 
remains at rest when the slider i rearranges, and the re- 
lease of the elastic energy loaded in the spring leads to a 
change 

AFi^2 = -2 (-I)XMs,, - f^d)N (A4) 

in the force Fi^2 between the sliders -see Eq. ([7]). The 
value of ^s,i needs to be drawn anew from p(/is), the 
other frictional coefficient remaining unchanged, and the 
process is here again iterated, leading back to Eq. (jAip . 
as long as a macroscopic motion does not emerge. 

At the end of this iteration process, the critical value 
tanac has been obtained, and the evolution on the hor- 
izontal substrate resumes (we recall that the evolution 
on the tilted substrate is only a 'virtual' test, that does 
not influence the 'real' horizontal evolution). Obviously, 
the obtained value of tanac depends on the values of 
{Fi^2, Ms,i, A^s,2, 0} given at the beginning of the test, 
which are all stochastic variables. It is thus necessary to 
average tanac over the dynamics. Since, on the other 
hand, we wish to obtain the time-dependence of the crit- 
ical angle to study aging effects, we proceed through a 
time scale separation as follows: tanac is averaged over 
all the rearrangements occurring in a given time step of 
duration Tth, and also over a large number of indepen- 
dent realizations of the horizontal dynamics, with ran- 
dom initial conditions. Then the time-dependence on a 
scale larger than Tth remains, and the time-step averages 
(. . .)t^ are labeled by the corresponding time tg = qrth, 
with q integer, yielding the effective static frictional co- 
efficient ^J,s^{tq) = (tanac)t, . 
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